Simulations of the 3d Ginzburg-Landau Model with Soft Amplitudes 
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Using cluster Monte Carlo simulations, the 3d complex Ginzburg-Landau model reveals a first 
order transition when the amplitude \tjj\ of the complex field ip is sufficiently soft, i.e. adapts itself 
to the phase configurations of the field. This transition is driven by phase fluctuations in agreement 
with a previous analytical approach. 
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The 3d complex Ginzburg-Landau (GL) model (for a 
review see jjjjpj) originally introduced by Landau to de- 
scribe general features of phase transitions, is still of great 
relevance in domains like superconductivity, liquid crys- 
tals or particle physics, for example. In particular, the 
role of amplitude fluctuations and the order of the tran- 
sition is still a subject of debate especially when one al- 
lows large amplitude fluctuations. The £-expansion || in 
the framework of the renormalisation group predicts that 
the 3d GL model should belong to the same universality 
class as the 3d XY model. On the other hand, a recent 
variational approach [Q has revealed the possibility of a 
first order transition induced by phase fluctuations. In 
this scenario, the interplay between phase and amplitude 
leads the amplitude to adapt itself to the configuration of 
phases and vice-versa. When the amplitude is soft, then 
the adaptation is so strong that a first order transition 
replaces the usual continuous transition. 

When two variables are coupled, it is well-known that 
fluctuations can change the order of the transition. For 
example, fluctuations of the magnetic field change the 
GL transition to a first order transition for type I su- 
perconductors ||. The three state Potts model in two 
dimensions is an opposite example: mean field theory 
predicts a first order transition, whereas the actual tran- 
sition is continuous. For the GL model, it has been pre- 
dicted |i] that the transition could be first order if the 
amplitude of the field ip — |-0|e 1 ^ is sufficiently soft. If 
this is the case, phase fluctuations would drive the tran- 
sition via the amplitude to a first order transition. It has 
been shown Q that amplitudes are dominated by phase 
fluctuations in the domain of the first order transition 
for d = 3. Therefore amplitude fluctuations can be ne- 
glected in this regime and phase fluctuations drive the 
XY transition to a first order transition via their action 
on the amplitude. 

Bormann and Beck Q have also shown that amplitude 
fluctuations might alter the cooperative phenomenon oc- 
curring with phases, in particular in dimension 2. Like 
the XY model, corresponding to a fixed value of the am- 
plitude, the 2dGL model can be mapped onto a Coulomb 
gas describing vortex-antivortex pairs. As soon as one 
allows for amplitude variations, these topological excita- 
tions become energetically more favorable. Taking into 
account gaussian amplitude fluctuations, Bormann and 
Beck have shown that the system may be driven into a 



regime where - according to Minnhagen's phase diagram 
- a first order transition replaces the usual Kosterlitz- 
Thouless scenario. 

Monte Carlo simulations have become increasingly im- 
portant to study statistical systems. Great progress has 
been made by using non-local algorithms for spin sys- 
tems. Swendsen and Wang [|| have used the Fortuin- 
Kastelyn || percolation mapping for the Potts model 
to define domains or clusters to be inverted with zero 
free energy cost. Cluster Monte Carlo simulations have 
a small critical slowing down comparing to standard 
Metropolis algorithms and therefore allow to study the 
critical region for large systems. 

Concerning the 3d GL model, to our knowledge, there 
is no systematic study of the influence of amplitude fluc- 
tuations. Most of the existing simulations (for example 
JlOt) are made in a domain where the XY behaviour dom- 
inate the transition. Here we have combined a standard 
Monte Carlo algorithm for the amplitude and a Wolff 
algorithm for the phase <fi in the same spirit as for 
the $ 4 model where the real field <E> is the product 
of its amplitude and of a discrete variable s = ±1. The 
cluster Monte Carlo algorithm allows to reduce the sim- 
ulation time in the region where the transition is XY-like 
or weakly first order. Indeed in this regime, the corre- 
lation length is very large, and therefore, using a clus- 
ter algorithm reduces strongly the critical slowing down 
ill- 

The aim of this letter is to verify the analytical pre- 
diction Q of the first order nature of the GL transition 
when amplitudes are soft by using cluster Monte Carlo 
simulations. 

According to Ginzburg-Landau theory, we define the 
effective hamiltonian functional 



H[ip] = J d d r 
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where a, b and 7 are coefficients derived from a micro- 
scopic model. We take a < so that the critical region 
can be reached by varying the temperature in the Boltz- 
mann factor. We put our system on a lattice, with lat- 
tice spacing e. In order to establish the phase diagram, 
we normalize the hamiltonian by setting -ip = %p / yj\a\/b, 
u = r/£, where £ 2 = 7/|a| is the mean field correlation 
length. The normalized hamiltonian is then: 
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where we have removed the constant term, and /i points 
to the <i nearest neighbours of the site i. N is the total 
number of sites. Only two competing parameters remain: 



a := e*/e 
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ct controls amplitude fluctuations. When cr is large, then 
the amplitude \\p\ is forced to take the value 1 to avoid too 
large costs for the energy and the only physics is then the 
one of the XY model. Vq corresponds to the zero temper- 
ature phase stiffness and is proportional to the superfluid 
density |ctj/6 in the language of superconductivity. When 
Vb/Tc is large, the critical region is small and the system 
has the properties of the mean field theory. When Vq /T c 
is of the order of 1 , phase fluctuations become very large 
and yield a critical temperature which is an upper 
bound for the true critical temperature T c fl3f| . 

The canonical partition function Z is then the trace 
over all possible configurations of the complex field tp: 



Z = Dip e 



-H/T 



(3) 



where the reduced temperature is defined as T = T /Vq. 

In the Swendsen-Wang algorithm, many clusters are 
generated and then flipped. Here, our idea is to use a 
standard Monte Carlo procedure for the amplitude 1^1 
and a Wolff algorithm for the phase cp as it has been 
done for the real <i> 4 model jjjj . A Wolff algorithm is sim- 
ilar to the Swendsen-Wang algorithm but only one single 
cluster is flipped at each step. The detailed procedure is 
the following: 

1. With fixed amplitude \tp\, introduce phase variables 
4> of the complex field ip, and take a site k (the 
seed) at random. Choose a random direction z with 
z| = 1. For a site i, we define a unit vector spin 
Si = {Re[tpi},lm[tpi]}/\tpi\. Nearest neighbours Sj 
of the cluster can be added to the cluster only if 
the component Sj ■ z has the same sign as the seed 
component Sfc • z (i.e. they point to the same half 
plane perpendicular to z). Then, if Sj is already in 
the cluster and a nearest neighbour Sj of s, points 
the same half plane, then Sj is added to the cluster 
with probability 



p = 1 - exp [2(z • Sj)(z • Sj)] 



(4) 



If the cluster stops growing, the process is inter- 
rupted and all spins of the cluster are flipped with 
respect to the axis perpendicular to z: 



2(z 



(5) 



Considering two configurations fi and v with an en- 
ergy difference AH, the ratio of the selection prob- 
abilities W(fi — > v) and W(v — > fj,) is then equal 
to e~$^ R , and the detailed balance is fulfilled (see 
the article of U. Wolff for the proof Jill). 

2. Adapt the amplitude \ip\ to the phase configuration 
with a standard Metropolis Monte Carlo algorithm. 
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FIG. 1. Comparison between standard MC simulations 
(lines) and cluster MC simulations (points). (a = 18(T — 1), 
6= 18) 

We present now the results of the cluster Monte Carlo 
simulations that have been performed on the isotropic 3d 
GL model with a complex field tp. The expectation value 
of a variable A is defined by: 



(A) 
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where M is the number of MC sweeps, and A m the value 
of A at the end of the sweep m. One sweep consists of 
points 1) and 2) explained above. 

The first measured quantity is the mean square ampli- 
tude: 
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In order to determine the critical temperature, the helic- 
ity modulus is also computed: 



r '' = \ Y1 \^i\\^i+^\cos((j>i - <j> i+ll ) 



1+jLl I 



(8) 



The helicity modulus is a measure of the correlation 
of phases in the direction /i. is zero above T c and 
different from zero below T r . 
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FIG. 2. 
values of a. 
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Points are for TV = 10 and dashed lines for 



N — 15 . The slanting dashed line shows the XY transition 
ending at T Ct xY ~ 2.2 (cr — > oo). 

A quantity of interest is also the XY energy of the 
phases. As in the XY model it is: 




(9) 



t )]/3 is the normalised 



where fa 
XY energy. 

Using different lattice sizes, enough sweeps have been 
used in order to have error bars of the order of the sym- 
bols used to plot the curves. In practice, the thermal 
equilibration was reached between 200 and 1000 sweeps, 
and then 10 4 sweeps were sufficient to have the required 
precision for N = 10 3 . We used lattice sizes from 3 3 to 
15 3 which is enough to reveal the first order transition. 
For point 2), we let \ip\ fluctuate from to a limit L 
that is chosen so that the probability of the least proba- 
ble configuration is less than 10~ 10 . This choice gives an 
acceptance rate of new amplitude configurations that is 
approximately between 30 % and 50 %. 

In order to control the validity of this method of cluster 
Monte Carlo simulations, we compare it with standard 
Monte Carlo simulations performed on the same model 
by Nguyen and Sudb0 jl0| . In figure [j], one can see that 
although the system sizes arc different, 15 3 for the cluster 
MC method and 60 3 for the other one, the two methods 
give results that are quite close. 

Results for different a are shown in figure || for N = 10 3 
(points) and N = 15 3 (dashed lines). The mean square 
amplitude (IV'I 2 ) is shown as a function of the reduced 
temperature T. The intersection between (l^) 2 ) and the 
slanting dashed line shows the XY critical temperature 
as it would be obtained with the XY model and a non- 
fluctuating coupling constant K(f) = (|^| 2 }/T. At T c , 
the we have the equation: 



In order to control the limiting case of the XY model, 
simulations have been made for a = 100: the critical 
temperature T c (where = 0) is near the XY critical 
temperature T c xy ~ 2.2, and moreover T c obeys equa- 
tion (|i~(i|). The most striking feature is the appearance 
of a first order transition below a critical value a* rj 2. 
If a > cr*, there is no observable first order transition 
and the transition is completely driven by the phase as 
also shown by Nguyen and Sudb0 jlo|. When a < cr*, 
then < \ip\ 2 > exhibits a jump at the critical tempera- 
ture T c . Simulations made for larger systems also con- 
firm the presence of this critical point a*. Indeed, re- 
sults for N = 15 3 are very close to the curves obtained 
for N = 10 3 . We have also used standard Metropolis 
Monte Carlo simulations for the regime where the transi- 
tion is strongly first-order (cr << a*). The same results 
are obtained when comparing to the cluster Monte Carlo 
method although the Metropolis algorithm takes much 
longer simulation time than the cluster algorithm even 
in the first order regime. 

The thermal equilibration is faster for systems with 
small a. Therefore, we can exclude the possibility that 
the Monte Carlo trajectory in the phase space jumps to 
another value if one would take much more MC sweeps. 



2000 



1500 



1000 



500 









E \ 



-0 . 2 



-0 . 1 



. 1 



0.2 



. 3 



K(f c ) = (\j\ 2 )(f c )/f c = f c , XY . 



(10) 



FIG. 3. Energy probability distribution near T c as a func- 
tion of the energy E. (a = 1.5, f = 0.4945, TV = 4 3 ) 

A typical plot of the energy distribution P(E) near 
T c is shown in figure [| when the transition is first-order 
(cr = 1.5). The two peaks of the distribution have the 
same height, and this reveals two coexisting phases near 
T c . This plot was obtained by sampling 50000 sweeps 
for N — 4 3 . Larger system sizes only emphasize the two 
peaks. However, it is difficult to control the shape of 
the peaks because the first order transition is not due a 
fundamental deg eneracy in the hamiltonian like for the 
Potts model |fyj| . These two peaks disappear when the 
transision is continuous unlike the Potts model so that it 
is not possible to determine the critical point a* by using 
the sampling method of Lee and Kosterlitz 14|. 

An interesting feature is the anti-correlation between 
(\'tp\ 2 ) and (/) during the MC simulation (see figure ||). 
This anti-correlation is stronger at T c and is due to the 
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FIG. 4. Monte Carlo trajectories of |-i/i| 2 (thick line) and /, 
(dashed line) showing the anti-correlation between these two 
variables (the parameters are the same as in figure |^) . 

coupling between amplitudes and phases. 

In figure ^, the value of the mean amplitude jump 
A(|^| 2 ) is plotted as a function of the parameter a. The 
analytical result from reference [Q is also shown (contin- 
uous line). Our calculations show that by increasing the 
size of the system the simulations get closer to the ana- 
lytical result. However, a precise value of a* is difficult 
to determine because of the finite size effects that round 
the jump at the transition. Nevertheless, we can say that 
<7* is restricted to the domain: 



2 < 



(11) 



In conclusion, cluster Monte Carlo simulations per- 
formed on the 3d complex Ginzburg-Landau model ex- 
hibit a first order transition for soft amplitudes when 
the parameter a is approximately smaller than 2. These 
conclusions are in good qualitative agreement with a pre- 
vious analytical work that predicts a first order tran- 
sition for a < 4.5. Therefore, if we take into account the 
rounding of the jump at the transition in the simulations, 
the critical point a* , that separate first order and contin- 
uous domains, should lie between 2 and 4.5. In order to 
have a better accuracy and a precise value of the critical 
<7*, one should use much larger systems and finite size 
scaling. 

Concerning a possible observation of this mechanism, 
the first order transition between nematic and smectic A 
seems to be a good candidate to measure the influence 
of phase fluctuations on the smectic A density \tp\. It is 
believed that the first order nature of this transition is 
due to the HLM mechanism [^J, but recent experiments 
fl5| show some differences between HLM predictions and 
measurements. It is therefore not clear which mechanism 
is responsible for the first oder nature of the transition. 
Same conclusions apply to the case of type I supercon- 
ductors where however a latent heat has never been mea- 
sured. 

As already mentioned in Q, high temperature super- 
conductors near the overdoped regime appear to be good 
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FIG. 5. Size of the amplitude discontinuity as a function 
of a. We compare the simulations (points) with the analytical 
result of reference 14 . 



candidates for an observation of a first order transition, 
as for example the cuprate b^Si^CaC^Cv 
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